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Abstract 

Quantum entanglement entropy has a geometric character. This is illustrated by 
the interpretation of Rindler space or black hole entropy as entanglement entropy. 
In general, one can define a "geometric entropy", associated with an event horizon 
as a boundary that concentrates a large number of quantum states. This allows 
one to connect with the "density matrix renormalization group" and to unveil its 
connection with the theory of quantum information. This renormalization group has 
been introduced in condensed matter physics in a heuristic manner, but it can be 
conceived as a method of compression of quantum information in the presence of a 
horizon. We propose generalizations to problems of interest in cosmology. 

1 Entanglement entropy in some relevant geometries 

Entanglement or nonseparability refers to the existence of quantum correlations between 
two sets of degrees of freedom of a physical system that can be considered as subsystems. 
It is natural that two (sub) systems in interaction are entangled an that they still are 
entangled after their interaction has ceased. Particularly interesting situations arise when 
two entangled systems become causally disconnected because of an event horizon. 



1.1 Introduction: entanglement entropy 

The entanglement of two parts of a quantum system can be measured by the von Neumann 
entropy. This is defined in terms of the density matrices of either part. Let us consider, for 
later convenience, one part as "left" or "interior" and another as "right" or "exterior", in 
a yet imprecise sense. Then, let us represent states belonging to the left (or interior) with 
small letters and states belonging to the right (or exterior) with capital letters. A basis 
for the global states (left plus right) is {\a)} ® By representing the ground state in 

this basis as 

|0> = X>aA|a>®|A), (1) 

aA 
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we define a coefficient matrix ip a A- Then we have two different density matrices: 

Trip^ip ' Tr ip*ip T 

Correspondingly, we have two von Neumann entropies: 

1. S R = -Tr A (p R \np R ) 

2. S L = -Tr a (p L \np L ) 

Now it is important to recall the "symmetry theorem" , which states that both entropies are 
equal, S R = Sl (this can be proved in several ways, the most popular one appealing to the 
Schmidt decomposition of the entangled state. The equality of entropies implies that they 
are associated with properties shared by both parts, that is, with (quantum) correlations. 

More generally, two non-interacting parts of a quantum system can be originally in 
respective mixed states. After their interaction, which we describe as an arbitrary unitary 
evolution of the composite system, the initial density matrix Pl®Pr has evolved to p' LR . It 
is easy to see that the partial traces p' L and p' R , in general, have von Neumann entropies S' L 
and S' R such that S' L + S' R > Sl + S R [1]. Of course, if the initial state is pure Sl = S R = 0. 



1.2 Field theory half-space density matrix 

Let us now consider the quantum system to be a chain of coupled oscillators. Moreover, we 
shall chiefly work in the continuum limit, where the concepts and mathematical expressions 
are more transparent, in spite of dealing with non-denumerable sets of degrees of freedom 
(we will return to a discrete chain to describe the density matrix renormalization group 
algorithm). The action for this model, namely, a one-dimensional scalar field, is 

A[<p(x, t)] = Jdtdx Q [{d t vf - (d x <p) 2 } - , (3) 

where tp is the field. 

Let us obtain a path integral representation for the density matrix on the half-line of 
a system that is in its ground state [2, 3, 4]. In the continuum limit, the half-line density 
matrix is a functional integral, 

p[<PR(x),<p' R (x)] = J D<p L (x)ifo[<p L (x),<p R (x)]i(>Z[<p L (x),<p' R (x)], (4) 

where the subscripts refer to the left or right position of the coordinates with respect to 
the boundary (the origin). Now, we must express the ground-state wave- function as a path 
integral, 

ij>o[<p L (x), ip R (x)\ = J Dtp(x, t) exp (-A[tp(x, t)\) , (5) 

where t G (— oo, 0] and with boundary conditions <p(x, 0) = (Pl(%) if x < 0, and <p(x, 0) = 
(f R (x) if x > 0. The conjugate wave function is given by the same path integral and 
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Figure 1: Geometric setting for computing p((p R , (p' R ) = (<p' R \ exp(— 2ir £)\ip R ), with angular 
Hamiltonian C and temperature l/(27r). 

boundary conditions but with t G [0, oo). Substituting into Eq. (4) and performing the 
integral over (Pl(x), one can express p(ip R ,ip' R ) as a path integral over ip(x,t), with t G 
(—00,00), and boundary conditions ip R (x,0+) = (p' R (x), ip R (x,0—) = <Pr(x). In other 
words, p{(fiR, <p' R ) is represented by a single path integral covering the entire plane with a 
cut along the positive semiaxis, where the boundary conditions are imposed. 

1.3 Angular quantization and Rindler space 

In Euclidean two-dimensional field theory, the generator of rotations in the (x, t) plane is 
given by 

£ = J dx(tT u - xT 00 ), (6) 

in terms of the components of the stress tensor computed from the action (3). To simplify, 
one can evaluate it at t — 0. Let us consider a free action (V = 0). In the Schrodinger 
representation, we should replace the momentum II = d t ip with H(x) = i5/5ip(x). How- 
ever, as in canonical quantization, one rather uses the second-quantization method, which 
diagonalizes the Hamiltonian by solving the classical equations of motion and quantizing 
the corresponding normal modes. Let us recall that, in canonical quantization, if we disre- 
gard anharmonic terms, the classical equations of motion in the continuum limit become 
the Klein-Gordon field equation, giving rise to the usual Fock space. Not surprisingly, the 
eigenvalue equation for C leads to the Klein-Gordon equation in polar coordinates in the 
(x, t) plane, 

The free field wave equation in polar coordinates, 

( A + + ^ + ™ I )*- = °. < 7 > 

can be solved by separating the angular variable: it becomes a Bessel differential equation 
in the r coordinate with complex solutions I±ie(mr), £ being the angular frequency. We 
have a continuous spectrum, which becomes discrete on introducing boundary conditions. 



3 



One of them must be set at a short distance from the origin, to act as an ultraviolet 
regulator [2, 3, 4], necessary in the continuum limit. 

Therefore, the second-quantized field is (on the positive semiaxis t = <^ = 0, x = r) 

r d£ b e I ie (mx) +b\l-. H (mx) 



where we have introduced annihilation and creation operators and where the term that 
appears in the denominator is just for normalization, to ensure that those operators satisfy 
canonical conmutations relations. There is an associated Fock space built by acting with 
b\ on the "vacuum state" . These states constitute the spectrum of eigenstates of £, which 
adopts the form C = J d££b\b( (where the integral is replaced with a sum for discrete £). 

The type of quantization just exposed was first introduced in the context of quantization 
in curved space, in particular, in Rindler space. Rindler space is just Minkowski space 
(therefore, not curved) in coordinates such that the time is the proper time of a set of 
accelerated observers. Its interesting feature is the appearance of an event horizon, which 
implies that the ground state (the Minkowski vacuum) is a mixed (thermal) state [5]. The 
connection with black hole entropy and Hawking radiation is explained in the next section. 

It is pertinent to note that the functions I±in(mx) have wave-lengths that increase with 
x. It is illustrative to represent a real "angular wave", 

KiAmx) = 7 — — — \Lp(mx) — I_iAmx)]. 

v ; 2 sinh(Trf) L V ; V n 

This solution is oscillatory for x < £/m, with a wavelength proportional to x, and decays 
exponentially for x > £/m (Fig. 1). The fact that the wave length vanishes at x = is to 
be expected from the Rindler space viewpoint, because it corresponds to the horizon. 

1.4 Black hole entropy 

Let us consider the Schwarzschild geometry in the Kruskal-Szekeres coordinates u, v, de- 
fined by 

uv = 16M 2 — -1 exp — -1 , - = exp 



,2M J 1 V2M J v 2M 
The horizon is given by u = or v = 0, like in the Rindler geometry of the previous section. 
Then, if we further define 

Z — u + v, T = u — v, 

these coordinates behave like the ordinary Minkowskian coordinates. Moreover, for small 
u or v (or the large M limit), the curvature can be neglected and the geometry becomes 
locally the one of Rindler space. Therefore, we could define radial coordinates Z = p cosh r, 
T = p sinh r, and perform a radial quantization like we did in the preceding section. 

Once established that the geometry near the black-hole horizon is locally the Rindler 
geometry of the preceding section, we can readily transfer the form of the density matrix 
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Figure 2: Real radial wave K 8i (x). 



of a scalar field therein, where we now ignore (trace over) the degrees of freedom inside the 
horizon. Hence, we can define a von Neumann entropy associated with this density matrix. 
Furthermore, in so doing, we can appreciate that the concept of black-hole entropy takes a 
new meaning: in addition to being of quantum origin, this entropy is related with shared 
properties between interior and exterior, namely, with the horizon. In addition, the radial 
vacuum is a thermal state with respect to the original Schwarzshild coordinates, giving rise 
to Hawking radiation [5]. 

2 Quantum information and RG transformations 

2.1 Information theory and maximum entropy principle 

The entropy concept appeared in Thermodynamics but only took a truly fundamental 
meaning with the advent of information theory. In this theory, entropy is just missing 
information, while information itself is often called negentropy. To recall basic definitions, 
the information attached to an event that occurs with probability p n is /„ = — lnp n . 
Hence, the average information (per event) of a source of events is S({p n }) = J2nPnI n = 

The previous definitions, given by Shannon in his theory of communication, seem un- 
related to thermodynamic entropy as a property of a physical system. However, according 
to the foundations of Statistical Mechanics on Probability Theory (the Gibbs concept of 
ensembles), a clear relation can be established, as done by Jaynes [6]. Jaynes made con- 
nection with the Bayesian philosophy of probability theory, in which the concept of "a 
priory" knowledge is crucial. Indeed, although the exact microscopic state of a system 
with many degrees of freedom may be unknown, one has some "a priory" knowledge given 
by the known macroscopic variables. This is a particular case of Jaynes' adaptation of the 
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Bayesian probability theory, wich postulates that the best probability distribution to be 
attributed to a stochastic event is such that it incorporates only the "a priory" knowledge 
about the event and nothing else. This postulate amounts to Jaynes' maximum entropy 
principle: given some constraints, one must find the maximum entropy probability distri- 
bution (density matrix, in the quantum case) compatible with those constraints, usually, 
by implementing them via Lagrange multipliers. In particular, more constraints mean less 
missing information and lead to less entropy. 

2.2 Quantum information 

The concepts of Shannon's classical theory of communication have quantum analogues 
[1, 7]. Nevertheless, the quantum theory of communication is richer (and less intuitive!). 
Indeed, the key new notion in the quantum theory is entanglement (already described in 
the foregoing): if a state (an event) is entangled with the environment, we have the type 
of purely quantum phenomena to which the EPR paradox is associated. 

Schumacher posed the problem of communication of an entangled state [7] . (The tech- 
nical name is transposition, since the copy of a quantum state is not possible: no-cloning 
theorem [1].) His conclusion is that the von Neumann entropy of the state is the quantity 
that determines the fidelity of the transposition: it is possible to transpose the state with 
near-perfect fidelity if the signal can carry at least that information. 

The fidelity is simply defined as the overlap between normalized states, namely, | ip') \ 2 . 
It is directly related with the natural distance in the space of rays, that is, the angle be- 
tween rays. The geometrical meaning of this distance is best perceived by considering the 
complex projective space of rays, where it is called the Fubini-Study metric [8]. 1 Schu- 
macher maximization of the fidelity uses the Schmidt decomposition of the entangled state 
[7]- 

2.3 Renormalization group and information theory 

The problem of transforming one quantum state into another while preserving its funda- 
mental features is reminiscent of an operation performed with the quantum renormalization 
group (in this connection, see Ref. [11]). This operation does not intend to reach near- 
perfect fidelity but it is desirable that it reach as much fidelity as possible. We shall see 
that a particular formulation of the renormalization group, namely, the density matrix 
renormalization group, comes close to applying the concepts of quantum information the- 
ory, and in a similar way to Rindler space quantization [13]. In fact, the construction of the 
density matrix renormalization group algorithm is based on the Schmidt decomposition, 
in parallel to Schumacher fidelity maximization. 

1 It is very interesting to realize that the distance between the distribution probabilities defined by the 
outcomes of possible measurements of observables that distinguish the states if) and tjj' coincides with the 
statistical distance in the classical sense of Fisher [9, 10]. 
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Let us consider an entangled quantum state \ip) for which we seek an optimal reduced 
representation (a sort of quantum coding). That is, we must find a projection 

= £ ^ |a> ® M - |$ = £ 4>iA \i) ® 1^4) (9) 

aA iA 

to a subspace spanned by reduced basis {i} such that the distance S — — \ip)\ 2 is 
minimized. This is the problem that S. White met in looking for an improved "real space" 
quantum renormalization group algorithm: In a renormalization group step, one must 
reduce the number of block states, but in a way that accounts for the influence of the rest 
of the system (that is, with no introduction of arbitrary boundary conditions) [12]. This 
amounts to consider the entanglement of the block and to formulate the above described 
problem. White found the solution in terms of the singular value decomposition, a well- 
known numerical algorithm, equivalent to the Schmidt decomposition of the entangled state 
and, hence, leading to discarding the smallest eigenvalues of the block density matrix; in 
other words, \i) are the eigenstates of p with largest eigenvalues. The proper formulation 
of the density matrix renormalization group in terms of quantum information concepts has 
been provided recently [14]. 



2.3.1 Density matrix renormalization group algorithm 

Method for a ID quantum system on a chain (e.g., a chain of oscillators): 

1. Select a sufficiently small "soluble" block [0,L]: I • • • 

2. Reflect the block on the origin: • • • I • • • 

3. Compute the ground state. 

4. Compute the density matrix of the block [0, L}. 

5. Discard eigenstates with smallest eigenvalues. 

6. Add one site next to the origin. 

7. Go to 2. 

One has to adjust this procedure in such a way that the iteration keeps the Hilbert space 
size approximately constant. The procedure can be performed algebraically for a chain of 
coupled harmonic oscillators [13]. Otherwise, it has to be performed numerically. In the 
continuum limit, the connection with the scalar field theory half-space density matrix 
previously calculated is clear. In particular, the property of "angular waves" of having 
vanishing wave-length at x = pointed out above provides an explanation of the form in 
which the density matrix renormalization group solves the boundary condition problem, 
that is, because of the concentration of quantum states near the boundary. To be precise, 
using the eigenfunctions of C instead of free waves, we have a basis in which the region 
close to x = (the boundary point) is more accurately represented than the region far 
from it when we cut off the higher i eigenfunctions. 
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3 Generalizations and applications 



3.1 Geometric entropy 

We have seen that the half line density matrix of a field theory has a geometric interpre- 
tation in Rindler space. Furthermore, the entropy of black holes can be understood as a 
generalization to a more complicated geometry. We may wonder if further generalizations 
are possible. 

In this connection, we recall the notion of "geometric entropy", introduced by C. Callan 
and F. Wilczek [4], as the entropy "associated with a pure [global] state and a geometrical 
region by forming the pure state density matrix, tracing over the field variables inside the 
region to create an 'impure' density matrix". Of course, their motivation was the earlier 
suggestion that back-hole entropy is of quantum-mechanical entanglement origin [2]. They 
proposed a generalization to different topologies but, actually, they only computed the 
Rindler space case, discussing the divergence of the entropy at the horizon [4] (the UV 
divergence of this type of entropy had been discussed in general in Ref. [2]). 

A different notion of geometric entropy can be deduced by purely geometrical means 
from the presence of horizons, namely, as associated with a spacetime topology that does 
not admit a trivial Hamiltonian foliation [15]. This type of topology prevents unitary 
evolution and leads to mixed states. 

In fact, it is only the second type of entropy that embodies the famous "one-quarter area 
law" for black holes, due to its origin in purely relativistic concepts (this was demonstrated 
for an earlier relativistic notion of entropy in Ref. [16]). On the contrary, the quantum 
notion of geometric entropy for a field theory involves UV divergences and needs renormal- 
ization before a comparison with the relativistic notion can be made (see the discussion in 
Ref. [17]). 

3.2 Application to cosmology 

The generalization of the concept of black-hole entropy to the de Sitter space and, hence, 
to cosmology is relatively old [18]. Of course, the concept of entropy in its traditional 
thermodynamical sense has been crucial in explaining the dynamics of inflation (now a 
standard paradigm), namely, in accounting for the reheating process (entropy generation). 
However, the relations between the traditional view and the one associated to quantum 
entanglement may lead to further insight, when they are properly formulated in the context 
of quantum information theory. 

In particular, we can regard the generation of entropy and fluctuations in de Sit- 
ter spacetime as a fundamentally quantum process leading to the celebrated Harrison- 
Zeldovich scale invariant spectrum of Gaussian fluctuations. If we consider the initial 
state for inflation as a pure quantum state (that some theory of quantum gravity will 
hopefully characterize some day), the de Sitter space horizon induces decoherence of the 
modes which cross it, irrespective of the actual inflaton dynamics and, therefore, of the 
details of the reheating process. This decoherence consists of a randomization of the phases 



8 



of the present quantum fields and naturally produces thermal Gaussian fluctuations. More- 
over, the symmetry of de Sitter space implies that each mode has the same physical size 
as it crosses the horizon, leading to the Harrison-Zeldovich power spectrum: 



P(k) 
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oc k . 



Actually, this spectrum is given by the equipartition theorem at the corresponding Hawking 
temperature T = H/(2n) (H = A/3) [18]. 

In this reasoning, one could also consider the initial state to be mixed (e.g., a thermal 
state) and apply the argument for entropy growth exposed in subsection 1.1. 

One may wonder if a density-matrix type renormalization group could help with the 
dynamics. As long as the fluctuations are Gaussian, we are in a similar situation to the 
case of harmonic oscillators commented above, which makes any renormalization group 
superfluous. However, as is well known, the gravitational instability produces non-linear 
evolution and leads to phase correlations. Therefore, the ideas presented here may be 
useful in setting up a renormalization group for the study of non-Gaussian fluctuations 
and its non-linear evolution. In particular, the Wilson or exact renormalization group 
irreversibility properties and can be connected with other methods of analysis of the non- 
linear evolution (for preliminary steps in that direction, see Ref. [19]). 
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